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ABSTRACT 

Cosmological hydrodynamical simulations of galaxy formation in representative re- 
gions of the Universe typically need to resort to subresolution models to follow some 
of the feedback processes crucial for galaxy formation. Here, we show that an energy- 
driven outflow model in which the wind velocity decreases and the wind mass loading 
increases in low-mass galaxies, as suggested by observations, can produce a good match 
to the low-mass end of the observed galaxy stellar mass function. The high-mass end 
can be recovered simultaneously if feedback from active galactic nuclei (AGN) and a 
correction for diffuse stellar light plausibly missed in observations are included. At the 
same time, our model is in good agreement with the stellar mass functions at redshifts 
z = 1 and z = 2, and with the observed redshift evolution of the cosmic star forma- 
tion rate density. In addition, it accurately reproduces the observed gas to stellar mass 
ratios and specific star formation rates of galaxies as a function of their stellar mass. 
This agreement with a diverse set of data marks significant progress in hydrodynam- 
ically modelling the formation of a representative galaxy population. It also suggests 
that the mass flux in real galactic winds should strongly increase towards low-mass 
galaxies. Without this assumption, an overproduction of galaxies at the faint-end of 
the galaxy luminosity function seems inevitable in our models. 
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1 INTRODUCTION 

The formation and evolution of galaxies is regulated by 
a complex interplay of accretion, star formation, and 
associated feedback processes. Unravelling the detailed 
physics involved is one of the greatest challenges in mod- 
ern astrophysics. High-resolution cosmological simulations 
guided by increasingly accurate observations have be- 
come one of t he most important tools for tackling this 
problem (e.g. Ijenkins et ail Il998l; ISpringel et all l2005bl ; 
iBovlan-Kolchin et al.l 120091 : Klypin et all l201ll ). However, 
even with today's fastest supercomputers, it is numerically 
still infeasible to fully resolve all the relevant processes that 
govern the formation of stars in galaxies and their back- 
reaction on the interstellar medium (ISM), even though 
remarkable progress towards this goal has recently been 
made for simulations of isola ted galaxies (|Hopkins et al.l 
l201ll . l2012l : lGuedes et al.ll201ll ). Thus, hy drodynamical sim- 



ulations of full gal a xy populations (e.g . Crain et "all 120091 : 



ISchave et ail |2010| ; IVogelsberger et al.r[20l3 r usually 



ploy simplified sub-resolution models to capture star for- 
mation and associated feedback processes in an approxi- 
mate way, with parameterizations that are guided by ob- 
servations. So-called semi-analytical galaxy formation mod- 



els (e.g. IGuo et all 1201 ll ) go yet a step further and also re- 
place the hydrodynamics itself with a highly approximate 
schematic treatment, using as input halo merger trees ob- 
tained from coUisionless simulations. 

An important difficulty that every theoretical model of 
galaxy formation faces within the ACDM cosmogony is to 
explain the very different shape of the dark matter halo mass 
function and the galaxy stellar mass function (GSMF). The 
halo mass function exhibits an almost featureless power-law 
with a steep logarithmic slope of d log n /rilogM ~ —0.9 fo r 
halo masses below ~ 10 13 /i _1 M Q (e.g., iTinker et al.ll2008h . 
The GSMF, instead, shows a relatively flat slope at the low 
mass end, dlog (j>/d\og M ~ —0.2, and a break (known as the 
'knee') at a characteristic stellar mass of ~ 5 x 10 10 /i _1 Af Q , 
with an exponential suppression of t he abundance of still 
larger galaxies (e.g., lLi fe White! l2009h . 

Cosmological hydrodynamical simulations that account 
for radiative cooling but lack feedback processes associ- 
ated with star formation fail to reproduce this behaviour 
and strongly overpred ict the stellar mass fractions of ha- 
los (iBalogh et alj|200ll ). It has hence been realized early on 
that some kind of feed back mechanism needs to be included 
l|White fc Frenklll99ll ). However, as the relevant physical 
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scales of the multiphase ISM or of AGN can not be re- 
solved in large volume simulations, previous work has intro- 
duced a number of phenomenologica l sub-resolution mod- 
els fo r treating such feedback effects JSpringel fc Hemauistl 
20031; ISchave fc Dalla Vecchial l200Sl; ISpringel et all 12005a 



Borgani et all 120051 ; ISiiacki et al.l 120071 ; iBooth fc Schavel 



2009). In most works, the primary feedback physics invoked 



has been supernovae, which are thought to be crucial for 
reducing the cosmic star formation rate to reasonable lev- 
els. However, different approaches were used in numerical 
models to account for the feedback energy supplied by su- 
pernova explosions. Many implementations inject at least 
part of the energy thermally into the ISM to counteract ra- 
diative cooling losses. However, in its simplest form, such 
an approach is relatively ineffective and does not create 
outflo ws, prompting a number of 'delayed cooling' models 
(e.g. iThacker fc Couchmanl [2000; Sca nnapieco et all I2006T ) 
that try to suppress a rapid cooling loss of the injected 
heat energy. Alternatively, a kinetic feedback component 
which directly drives galact ic winds has often been invoked 
i|Springel fc Hernquist|[2003l ). 

In the simplest parameterization, winds are driven with 
a fixed mass loading r), i.e. by assuming a fixed ratio between 
wind mass flux and the galaxy's star formation rate. If the 
assumed amount of energy avai lable for the kinetic feedback 
per unit mass of stars formed (|Springel fc Hernquistl [20031 ) 
is also kept fixed, a constant wind velocity results. Hence 
we will call such models constant wind models. They were 
found to significantly reduce the abundance of galaxies with 
stellar masses below a few times 10 10 Mq, but the strength 
and exact mass-scale of the suppression d epends on the de- 
tailed choice of the wind parameters (e.g. ICrain et al.ll2009l . 
see also our Fig. [T}. These simulations have, however, so far 
neither been able to reproduce the slope of the low-mass end 
of the observed GSMF nor its overall shape and normaliza- 
tion. Furthermore, simulations that did not include AGN 
feedback tended to overpredict the abundance of massive 
galaxies. 

In comparison, semi-analytical models have been much 
more successful in reproducing the global properties of 
galaxies, including their abundance as a function of stellar 
mass (i.e. the GSMF) and the redshift evoluti on of primary 
characteristics of the ga laxy population (e.g.. ICroton et all 
l200d ; lB~ower et al.l2006t ). From such studies it was concluded 
that the faint-end slope of the galaxy luminosity function 
can be made compatible with observations if the ratio of 
the mass flux of gas ejected from galactic discs by super- 
nova and stellar winds to the star formation rate is scaled 
with properties of the host galaxy. This has often been pa- 
rameterized by scaling the ejection efficiency with an inverse 
power of the disc rotation vel ocity, ~ v 7" , with typical val- 
ues for a in the range 2 t o 3.2 l|Cole et al. 2000; B enson et al.l 
120031 ; iBower et al.ll2012T ) . It was also found that by including 
simple models for AGN feedback the abundance of massive 
galaxies could be significa ntly reduced and made compati- 
ble w ith observations (e.g.. ICroton et a l. 2006; Bower et alj 
l2006h . 

Observationally, the mass loadings and velocities of 
galactic winds are still not accurately constrained, but some 
progress has been made in recent years by measuring Na , 
Na 1 and Mg 11 absorption lines in outflows. iMartinl (|2005l ) 
discovered a correlation between the wind velocities and star 



formation rates of galaxies, v w ~ SFR 0,35 . For galaxies with 
measured rotation curves, a positive correlation between the 
galaxies' circular velocities and the terminal wind velocities 
was found, and maximum wind velocities slightly smaller 
than the galaxies' escape velocities were inferred. IMartinl 
(2006) estimated mass outflow rates in the winds' cold com- 
ponents of the order of 10% of the SFR. The achieved level 
of accuracy did, however, not allow firm conc lusions about 
how th e mass loading scales with wind velocity. I Weiner et al.l 
(2009) observed a similar correlation between SFR and out- 
flow velocity, v w ~ SFR ' 3 , in z ~ 1.4 galaxies, as well 
as maximum win d velocities close to the escape velocity. 
IChen et al.l (|2010T ) found weak (yet noisy) positive correla- 
tions between outflow velocity and star formation rate den- 
sity or galaxy stellar mass, respectively. 

Motivated by these observations and the successes of 
semi-analytical models, first attempts have been made to 
include phenomenological wind models in cosmological hy- 
drodynamical simulations in which the wind velocity and 
mass loading scale with galaxy prop erties. We will refer to 
such m odels as variable wind models. lOppenheimer fc Pavel 
(2006) introduced a scaling of wind properties with galaxy 
velocity dispersion a by using the gravitational potential 
(f> at the location of a galaxy as a proxy for a, based on 
the approximation o ~ \/—l /24>. Assuming tha t winds are 
momentum driven (see, e.g.. lMurrav et aTll2005l ). they then 
adjusted the wind velocity v w and wind mass-loading r\ in 
their simulations according to w w ~ a and r\ ~ l/w w ~ 1/er. 
Compared to previous wind models this resulted in more 
realistic metallicities of the IGM as constrained by C IV ab- 
sorption studies. In addition, the observed cosmic star for- 
mation rate density at hig h redshift is reproduced better 
(see also Schav e~et al.ll2010T ). 

lOkamoto et al. I |2010T ) assumed an energy-driven wind 
scenario in which the wind properties scale according to 
u w ~(j and j) ~ l/v% ~ 1/<t 2 . Their model predictions were 
shown to be in broad agreement with the observed lumi - 
nosity function of Milky Way satellit es. iPave et al.l (l2011bT) 
used an updated version of the lOppenheimer fc Pavel 
(2006) momentum-driven model in which the galaxy veloc- 
ity dispersion is instead estimated based on halo masses 
computed by an on-the-fly friend-of-friend group finder 
l|Oppenheimer fc Pavell2008l ). They demonstrated that this 
model c an roughly reproduce the faint-end of the GSMF 
(also see lOppenheimer et alJlsiOlCl : IChoi fc Nagaminell201lh 
and the specific star formation rates of galaxies in the 
10 10 - 10 11 Mq stellar mass range. Furthermore, their vari- 
able wind simulations match the slope of the gas-phase 
metallicity - stellar mass relation better, and produce more 
realistic gas mass to stellar mass ratios in massive galaxies 
jPave et al.ll2011ah . 

Overall, variable winds are clearly more successful in 
suppressing star formation in low mass galaxies due to the 
boosted wind mass loading in these objects, as well as in 
predicting more realistic ISM masses and metallicities. In 
this work, we will hence explore variable wind models fur- 
ther. We will outline our numerical methods and describe 
our simulation set in Section^ discuss how the z — GSMF 
depends on the details of the wind model in Section T3. II and 
investigate the impact of AGN feedback and the influence 
of diffuse light on the abundance of massive galaxies in Sec- 
tion 13.21 In Section 13.31 we assess whether the wind model 
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that reproduces the 2 = GSMF best can also match the 
cosmic star formation history, as well as the GSMFs at red- 
shifts z — 1 and 2. We finally explore the impact of winds on 
the gas masses and specific star formation rates of galaxies 
in Section [3. 41 and summarize our findings in Section [4] 



2 METHODOLOGY 
2.1 The simulations 

Throughout this paper, we make use of a set of cosmo- 
logical hydrodynamical simulations that were performed 
with the TreePM-SPH simulation code P-Gadget-3, an 
updated and s ignificantly extended version of Gadget-2 
|Springej|2o"o!h . T he simulations adopt the best-fit WMAP 
7-year cosmology l|Komatsu et al.ll201l[ ) with !Jm = 0.272, 
n A = 0.728, Q B = 0.0465, h = 0.704, and a 8 = 0.809. 

In this work, we want to explore the simulated galaxy 
populations and galaxy stellar mass functions both at the 
low-mass and the high-mass end. This calls simultaneously 
for a high mass resolution (to resolve the small objects) 
and a large volume (to find the rare high mass objects). An 
immense computational effort would be necessary to meet 
these requirements in a single simulation. To avoid this prob- 
lem, we have chosen to perform a set of simulations that 
cover a range of box sizes from 15 to 60/i _1 Mpc and a fac- 
tor of 64 in mass resolution. The parameters of all runs that 
we use in this work are summarized in Table [T] 

Strictly speaking, it is problematic to run simulations 
with a comoving box size of only 15 h~ Mpc until redshift 
z = 0. At low redshift the fundamental mode in such small 
simulation boxes becomes non-linear, but since a coupling 
to larger modes is absent, the evolution of these modes will 
not be fully correct. The simulation box may then not be a 
representative sample of the universe any more. We checked, 
however, that the quantities in which we are most interested 
in, i.e. the GSMF and the stellar mass fraction of halos, are 
not significantly affected by finite box-size effects other than 
small- number statistics at the high-mass end (compare, e.g., 
the solid green curves in Figs. and [S]). The 15/i _1 Mpc 
boxes can thus be reliably used to investigate the effects of 
winds on the low-mass end. For the high-mass end, we use 
simulations with a box size of 60 h~ 1 Mpc. 



2.2 The wind model 

Inspired by the observations discussed in Section[T] we want 
to study how galaxy formation is affected by galactic winds 
whose velocities and mass loadings scale with the properties 
of the galaxies form which they are launched. More specif- 
ically, in our simulations we assume that the wind velocity 
u» is directly proportional to a galaxies' escape velocity w csc , 
i.e. 

V„ = K, X V csc . (1) 

The choice of the proportionality constant k will be dis- 
cussed in more detail in Section [3~T1 The escape velocity from 
the centre of a galaxy that is well described by a Navarro- 
Frenk- White (NFW) profile |Navarro et al.lll996l ) is given 



by 

^c = «200eXyC i+ 2 c) C __^_ , ( 2 ) 

where c = r s /?-2ooc is the halo's concentration, i.e. the ratio 
of the NFW-scale radius r s to the radius r200c at which the 
mean enclosed density is 200 times the critical density of the 
Universe. V2oq c = \J GA^ooc/^ooc is the circular velocity at 
»"200c, where G is the gravitational constant, and M200C is the 
mass enclosed within radius r200c In our simulations, we run 
an on-the-fly friends-of-friends (FoF) group finder at regu- 
lar intervals and store for each gas particle the mass of the 
FoF-group it is part of. Based on this mass, the host halo's 
escape velocity can be estimated using Eq. @ and the mass- 
concentrati on relation c = 4.6 7x (M2oo c /fl0 14 fe~ 1 M w 1)~ 11 
obtained bv lNeto et al.l (120071) fro m the Millennium N-body 
simulation (|Springel et alJl2005bl )R 

Once we have fixed the velocities of galactic winds, we 
need a prescription for the amount of gas that enters into 
such winds. A simple choice that was used in most previous 
sub-resolution models of wind formation is to link the power 
that is available for driving the wind to the star formation 
rate of the galaxy from which it is launched. In particular, 
the wind energy has often been assumed to be a specific 
fraction of the available sup ernova feedback energy (see e.g. 
ISpringel &: Hernquistl l2~00St ). The latter is proportional to 
the star formation rate, with a proportionality constant that 
depends on the assumed initial stellar mass function (IMF). 
Under these assumptions, the mass flux M w of the wind is 
related to the galaxy's star formation rate M* by 

1-2 

-M w v w = esN.wM*, (3) 

where £sn,w is the energy available to the wind per unit mass 
of long-lived stars formed. For such energy-driven winds, the 
mass-loading r\ = M w /M* scales with r\ ~ v^ 2 . Thus, winds 
with velocities given by Eq. JTJ) will have significantly larger 
mass-loadings in low-mass halos and are therefore expected 
to suppress star formation more strongly there. 

In reality, the amount of energy that is available for 
driving a galactic wind may not only depend on the galaxy's 
star formation rate. For example, the fraction of injected 
energy that is lost by radiative cooling could depend on 
halo mass. On the other hand, momentum that is injected 
by stellar winds, supernovae and AGN cannot be radiated 
away. Accordingly, if galactic winds are momentum-driven, 
it would seem more natural to assume a proportionality be- 
tween the winds' momentum flux (rather than its energy 
flux) and the star formation rate of the host galaxy. In this 
case, the mass-loading would scale like n ^ v' 1 . 

In this work, we will explore the effects of both a ~ 
v~ 2 and a ~ v" 1 dependence of the mass-loading on the 
wind velocity, where we assume in both cases that the wind 
velocity is related to the halo's escape velocity by Eq. UJ. 
For reference, we will also show results base d on the galactic 
wind model originally introduced by ISpringel fc Hernquistl 



1 As the wind velocity is not very sensitive to the mass- 
concentration relation, we ignore redshift ev olution effects and 
use the z = relation from lNeto et aL I 1120071 1 also at z > 0. 
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simulation model name 


box size 
[/i-iMpc] 






[h-i-Me,] 


e 

[ft-ikpc] 




?/ 


/uv 


IMF 


L15N128csf 


15 


2 




128 3 


2 x 


10 7 


4.5 


- 


- 


1 


Sal. 


L15N128csf,cw 


15 


2 




128 3 


2.0 x 


10 7 


4.5 


484 kms -1 


2 


1 


Sal. 


L15N128csf,cw,polar 


15 


2 




128 3 


2.0 x 


10 7 


4.5 


484 kms -1 


2 


1 


Sal. 


L15N128csf,mdvwl.25 


15 


2 




128 3 


2.0 x 


10 7 


4.5 


1.25 X Vesc 


~ "w 1 


1 


Sal. 


L15N128csf,edvwl.25 


15 


2 




128 3 


2 x 


in 7 


4.5 


1.25 X Vesc 


~ «w 2 


1 


Sal. 


L15N128csf,edvw0.6 


15 


2 


X 


128 3 


2.0 x 


10 7 


4.5 


1.25 X Vesc 


~ «w 2 


1 


Sal. 


L15N128csf,edvw0.6,uv,chab 


15 


2 


X 


128 3 


2.0 x 


10 7 


4.5 


0.6 x v csc 


~ «w 2 


3 


Chab. 


L15N256csf 


15 


2 


X 


256 3 


2.5 x 


10° 


2.25 






1 


Sal. 


L15N256csf,cw 


15 


2 


X 


256 3 


2.5 x 


10° 


2.25 


484 kms - ^ 


2 


1 


Sal. 


L15N256csf,edvw0.6,uv,chab 


15 


2 


X 


256 3 


2.5 x 


10 6 


2.25 


0.6 X v csc 




3 


Chab. 


L15N512csf,edvw0.6,uv,chab 


15 


2 


X 


512 3 


3.2 x 


10 5 


1.125 


0.6 X v csc 


~ «w 2 


3 


Chab. 


L60N512csf 


60 


2 


X 


512 3 


2.0 x 


10 7 


4.5 






1 


Sal. 


L60N512csf,cw 


60 


2 


X 


512 3 


2.0 x 


10 7 


4.5 


484 km s" 1 


2 


1 


Sal. 


L60N512csf,edvw0.6,uv,chab 


60 


2 


X 


512 3 


2.0 x 


10 7 


4.5 


0.6 X D csc 


~ «w 2 


3 


Chab. 


L60N512csfbh,edvw0.6,uv,chab 


60 


2 


X 


512 3 


2.0 x 


10 7 


4.5 


0.6 x v csc 


~ «w 2 


3 


Chab. 



Table 1. Summary of the parameters of the simulations used in this work. Given are the comoving box size, the number of simulation 
particles n par t in the initial conditions (half of them are gas and half of them are dark matter particles), the gas particle mass m gas , 
the comoving gravitational softening (Plummer equivalent), the wind velocity ti w , the wind mass loading rj, the UV heating boost factor 
/uv, and the assumed IMF (Salpeter or Chabrier). 



which assumes a constant wind mass-loading r\ and 
wind velocity v w . 

Unless mentioned otherwise, we assume our winds to 
be isotropic, i.e. gas particles entering the wind near the 
centre of a galaxy will receive a kick w w in a random di- 
rection. We checked, however, that adopting polar outflows 
does not alter our resul ts significantly (see, e.g., Fig . [T]) , in 
line with the findings of ISpringel fc Hernquistl (|2003h . Once 
particles enter the wind through stochastic selection from 
the star-forming phase, they are briefly decoupled from hy- 
drodynamical interactions until their SPH density falls be- 
low 10% of the threshold d ensity for star formation (as in 
ISpringel fc Hernquistl [20031 ). To protect against rare cases 
in which wind particles do not reach a region of low enough 
density (this could happen if, e.g., the imposed wind velocity 
is very low), the hydrodynamical forces are turned on again 
after an elapsed time equal to 2.5% of the Hubble time at 
the simulation redshift. 



2.3 The AGN feedback model 

In a subset of our simulations, we include a model for 
the growth of supermassive black holes (BHs) and associ- 
ated feedback processes. Here, we briefly summarize how 
these processes are fol l owed. Full details a re given elsewhere 
|Springel et al.ll2005al : ISiiacki et aLll2007l V 

We assume that low-mass "seed" BHs are produced 
sufficiently frequently such that every halo above a cer- 
tain threshold mass can be expected to contain at least 
one such BH at its centre. In the simulations, we use the 
on-the-fly friends-of-friends group finder that is also used 
for our variable wind model to put seed BHs with a mass 
WBH.secd = 1O 4 /i~ 1 M0 into halos when they exceed a mass 
of 5x lO lo h _1 M and do not contain any BH yet. The BHs 
are represented by collisionless sink particles and are allowed 
to grow by mergers with other BHs and by gas accretion, 
assuming a Bondi-Hoyle-Lyttleton prescription for the ac- 
cretion rate, subject to an imposed upper limit equal to 



the Eddington limit. Two BHs are merged when they fall 
within their local SPH smoothing lengths and have suffi- 
ciently _jimallj-ela^^ velocit^^ previous 
work (|Springel et al. I l2005al : ISiiacki et al.1l2007l ), we do not 
consider any relative motion of a BH with respect to the sur- 
rounding gas when estimating its accretion rate. This was 
done here as the code does not follow the motion of the 
BHs self-consistently. Instead, they are effectively pinned to 
the centres of their host halos by recentering them on the 
potential minimum at every time step. 

Motivated by growing theoretical and observational ev- 
idence that AG N feedback is compose d of two characteristic 
modes (see e.g. IChurazov et al.ll2005l ). we use two distinct 
feedback models d epending on the BH accretion rate (see 
ISiiacki et al1l2007h . For large accretion rates above 0.01 of 
the Eddington rate, the bulk of AGN heating is assumed 
to be in the form of radiatively efficient quasar activity 
with only a small fraction of the luminosity being thermally 
coupled to the surrounding medium. We adopt this ther- 
mal heating efficiency to be 0.5% of the rest mass-energy 
of the accreted gas, which reproduces the observed rela- 
tion between BH mass and bulge stellar velocity dispersion 
l|Di Matteo et al.ll2005h . For accretion rates below 0.01 of 
the Eddington rate, we assume that feedback is in the so- 
called "radio-mode", where AGN jets inflate hot, buoyantly 
rising bubbles. The mechanical feedback efficiency provided 
by the bubbles is taken to be 2% of the accreted rest mass- 
energy, which is in g ood agreement with X-ray observations 
of elliptical galaxies (|Allen et alj|2006h . ISiiacki et al.l (|2007h 
showed that this bubble model leads to a self-regulated resid- 
ual BH growth in large halos and yields realistic BH mass 
densities. It was also found that feedback from the BHs sig- 
nificantly lowers the baryon fractions and X-ray luminosities 
of poor galaxy clusters and groups, bringing them into ex - 
cellent agreement with observations l|Puchwein et al.| [2008). 
Furthermore, the luminosities of central group and cluster 
galaxies are reduced an d made compatible with observations 
i|Puchwein et al.ll2010h . 
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2.4 Computing galaxy stellar mass functions and 
accounting for diffuse light 

In our simulations, we find galaxies and their stellar masses 
by applying the Subfin d substructure finder (|Springel et al.l 
l200ll ; lDolag et al.ll2009l ) to FoF groups identified with a link- 
ing length of 0.2 times the mean interparticle spacing. Sub- 
find decomposes each FoF group into a main halo and self- 
bound substructures. We consider the stellar content of each 
of these objects to be a galaxy and compute its stellar mass 
by simply summing up the masses of the corresponding star 
particles. This allows as to compute GSMFs. 

In our analysis, we also examine stellar to halo mass 
ratios, based on halo central galaxies. In our simulations, 
the latter correspond to the star particles that are part of 
a FoF groups' main (sub)halo. The group's M200C value is 
adopted as halo mass and is based on a spherical overdensity 
mass definition^. 

Galaxy groups and clusters do not only contain stars 
that are bound to their member galaxies but also a 
signifi cant diffuse stellar co mponent (e.g., IZibetti et alj 



120051 ; iMcGee fc Balogbl |2010 | ) which is als o found in cos- 
molo gi cal simulations (|Murante et al.ll2004l; IWillman et alj 
20041; ISommer-Larsen et all 120051 ; iMurante et al.l 120071 ; 



Puchwein et al.ll2010l ). These diffuse stars are often not con- 
sidered to be part of a group's central galaxy, or are simply 
missed altogether in shallow observations. We will, there- 
fore, explore in Section \3. 21 how the massive end of our sim- 
ulated GSMF changes if we exclude such diffuse stars when 
computing the stellar masses of central galaxies. We do this 
based on a surface brightness threshold. More precisely, we 
compute the luminosities of the star particles in the simula- 
tions in the K and V-band based on their ages and metal- 
licities using the stellar population syn thesis model library 
GALAXEV (|Bruzual fc Charlotl 120031 ) . Next, we calculate 
the average surface brightness and enclosed stellar mass pro- 
files of galaxies in bins of halo mass. In this way, by comput- 
ing the profiles based on a large number of galaxies for each 
halo mass bin we avoid uncertainties in the surface bright- 
ness due to particle noise. We then derive for each halo mass 
bin the radius at which the profile falls below the adopted 
surface brightness threshold and the fraction of stellar mass 
within this radius. Finally, we correct the stellar masses of 
all galaxies in the corresponding halo mass bin by this frac- 
tion. However, in order to avoid reducing the stellar mass of 
low-mass galaxies which are not well resolved in the simula- 
tion, we do not remove stars that are less than one softening 
length from a galaxy's centre. 

We have checked that this procedure yields very similar 
results to computing SPH-like surface brightness estimates 
for all star particles, and then excluding those that fall below 
the threshold value and are at radii larger than one soften- 
ing length. For the threshold, we adopt values commonly 
used in the literature as a definiti on of diffuse light / BCG 
cutoff, i.e. V = 26.5 m ag/arcsec 2 jRudick et al.l 12000 ) and 
K = 20.7 mag/arcsec 2 (|Lin fc Mohr 2004)0 



2 -^200c is the mass within a spherical region around the halo 
centre in which the mean density is 200 times the critical density 
of the Universe. 

3 The BCG luminosity estimates in iLin k. Mo hr (2004) corre- 
spond to an effective cutoff at K ~ 20.7mag/arcsec 2 . 



3 RESULTS 

3.1 The low-mass end of the galaxy stellar mass 
function 

We first explore how the details of the wind model affect 
the GSMF and the stellar mass fractions of halos at red- 
shift z = 0. Figure [T] illustrates the effects of constant 
winds with commonly used settings for mass loading and 
wind velocity (77 = 2 and w w = 484kms _1 ). We com- 
pare runs wi th isotropic and polar cons tant winds (imple- 
mented as in lSpringel fc Hernquistll200o1 ) to a reference run 
that does not include any subresolution wind model. Ob- 
servational constra ints on the GSMF based on SDSS data 
|Li fc White 2009) and halo stellar fractions obtaine d by 
the abundance matching technique (|Guo et al.l l2010l ) are 
also shown for comparison. When winds are not included, 
the GSMF of well-resolved galaxies (i.e. with a stellar mass 
above ~ 5 x 10 9 h~ 2 Mq) is a power-law with roughly the 
same steep slope as the halo mass function over the corre- 
sponding range. Accordingly, the stellar to halo mass ratio 
varies only mildly. It decreases slowly with increasing halo 
mass, but the suppression of the abundance of galaxies with 
stellar masses below ~ 5 x 10 9 h~ 2 Mq is a resolution effect 
(see also Fig. [4j. 

In runs with winds, the stellar mass fractions of halos 
are reduced by roughly a factor of 4 to 5, somewhat less 
at the high-mass end, where the constant veloci ty winds be- 
come inefficient in expelling g £LS, £LS also found bv lCrain et al.l 
(2009). Otherwise, the suppression factor does not strongly 
depend on halo mass, as expected for winds that are faster 
than the escape velocity and have a constant mass loading. 
Correspondingly, the GSMF is shifted towards lower stellar 
masses, whereas its shape is only modified at the high-mass 
end. We note that the isotropic and polar wind models yield 
very similar results, thus the geometry of the outflows does 
not seem to be very important. Overall, we confirm that con- 
stant wind models are not able to reproduce the observed 
faint-end slope of the GSMF and fail to explain the different 
shapes of halo and stellar mass functions. 

We thus next try variable wind models in which the 
wind velocity is proportional to the escape velocity of the 
galaxy from whic h the wind is launched, as suggested by 
observations (e.g., [Martin 200^)- We consider both models 
with a momentum-driven and an energy-driven scaling of the 
mass loading and wind velocity. We choose the normalisation 
of the mass-loading such that r\ is equal to 2 if the wind 
velocity matches the constant wind model, i.e. 

77 = 2 x 484kms -1 /« w ( momentum-driven ), (4) 

energy-driven \ , . 
Salpeter IMF ) ' ^ ' 

For the energy-driven model, this means that the mass- 
loading is set by the available supernova en ergy, where th e 
latter is computed assuming a Salpeter IMF (Salpeter 1955). 
For the moment, we adopt the value n — 1.25 for the ratio 
of wind velocity to escape velocity The choice of K will be 
discussed in more detail later. 

As expected and shown in Fig. [21 a stronger scaling of 
mass-loading with galaxy size results in a shallower slope 
of the GSMF. In the energy-driven wind model, the stel- 
lar mass fractions of halos are no longer constant in the 
7 x 10 10 - 10 12 M Q halo mass range. Instead, the fractional 



r) = 2 x (484kms _1 ) 2 /i£ 
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Figure 1. Impact of galactic winds with constant velocity and mass loading on the low- mass end of the GSMF (left panel) and the 
stellar mass fractions of halos (right panel) at z = 0. Shown are results for a reference run without kinetic wind feedback (L15N128csf) 
and for simulations with isotropic (L15N128csf.cw) and polar (L15N128 csf,cw, polar) winds with constant wind velocity v w = 484km s 
and mass loading r? = 2. A fit to the observed GSMF jLi fc Whit 9*120091) is shown for reference. Also shown arc constraints on halo stellar 
fractions obtained using the abundance matching technique llGuo et*^Ul201fJ) , Above the resolution limit, both runs with constant winds 
and without winds strongly overpredict the low-mass end of the GSMF, as well as the stellar fractions in low mass halos. 



stellar mass content declines for lower halo masses. This 
trend is also seen in the observational constraints, although 
with a steeper slope. We conclude from this that the energy- 
driven wind model appears most promising for being able 
to reproduce the observed shallow slope of the GSMF's low- 
mass end. We thus concentrate on this model and explore 
in the following how the choice of n as well as assumptions 
about the available energy affect its predictions. 

Observations bv lMartinl (120051 ) (see their Fig. 7) found 
maximum wind velocities somewhat lower than the galax- 
ies' escape velocities n w ~ 0.6 x v CBC . Motivated by this re- 
sult we investigate whether slower energy-driven winds with 
k — 0.6 improve the agreement with observations. Figure [3] 
shows GSMFs and stellar mass fractions for this model as 
well as for the previously studied k = 1.25 simulation. The 
figure illustrates that slower winds with a correspondingly 
larger mass loading are more effective in suppressing the 
abundance of low-mass galaxies, while the number of galax- 
ies around the knee of the GSMF increases in such models. 
As a consequence, the « = 0.6 energy-driven model is able 
to reproduce the observed slope of the low-mass end of the 
GSMF. It even produces a knee at roughly the right position, 
however with a somewhat too shallow slope above the knee. 
Unfortunately, the normalisation of the predicted GSMF is 
still about a factor of 2 too large. 

An obvious way to reduce the overall normalisation 
of the GSMF is to increase the amount of available feed- 
back energy. In particular, we note that the amount of en- 
ergy injected by supernovae is not accurately known, leaving 
room for a potential upward revision of the available energy. 
Furthermore, our estimates thus far have been based on a 
Salpeter IMF, which implies a rather low fraction of stars 
that are massive enough to explode as a supernov a. Using a 
more recent calibration of the IMF like the one bv lChabrierl 
(2003) implies a factor of 1.8 larger number of supernovae 



per unit mass in stars formed. Adopting this revision, the 
mass-loading of the energy driven wind model is then given 
by 

i o r, -l\2 , 2 ( energy-driven \ , . 

„ = 1.8 X 2 X (484kms ) /Vm ^ Mp j . (6) 

Another heating source that is most likely underesti- 
mated in the simulations discussed above is photoheating as- 
sociated with h ydrogen and heliu m re ionization. This is im- 
pleme nted as in lKatz et alj (|l996l ) and lSpringel fc Hernquistl 
(2003) in our current simulations. However, during a sud- 
den reionization event, the simulation time steps in low 
density regions may be too large and the assumption of 
photoionisation equilibrium may be in accurate (see, e.g., 
IPeeples et al.l201 0; Puchwein et al . 2012, for a more detailed 
discussion), causing an underestimate of the energy injected 
during reionization. Furthermore, the photoheating rate is 
quite sensitive to the assumed spectrum of the cosmic ultra- 
violet background which is not well constrained. As a simple 
fix that allows us to get IGM temperatures that are com- 
patible with those inferred from analyses of the Lyman-a 
forest, we try boosting the photoheating rates in the code 
by a factor of 3. 

Both, assuming a Chabrier IMF and increasing the pho- 
toheating rate, lowers the normalisation of the GSMF and 
brings it into very good agreement with observations. This is 
illustrated in Fig. [3] The curves corresponding to simulations 
that assume a Chabrier IMF and a boosted UV heating are 
labelled Chab and UV in the legend of Fig. [3] respectively, 
as well as in all other figures of this paper. Checking the 
impact of these two changes individually, we found that in- 
creasing the supernova rate results in a larger reduction of 
the normalisation of the GSMF compared to boosting the 
photoheating rate. 

To summarize, we have now found a simulation model 
that accurately reproduces the low-mass end of the GSMF 
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Figure 2. Impact of scaling the velocity and mass loading of galactic winds with the host halo's escape velocity on the low-mass end of 
the GSMF (left panel) and the stellar mass fractions of halos (right panel) at z = 0. Shown are results for a reference run without kinetic 
wind feedback (L15N128csf), for simulations with isotropic winds with constant wind velocity d w = 484 kms -1 and mass loading n = 2 
(L15N128csf,cw), as well as for two variable wind models in which the wind velocity and mass loading scale with the host halo's escape 
velocity. In the latter models, the wind velocity is given by v w = 1.25 X v csc . The mass loading scales with the wind velocity according 
to n ~ v^j 1 for the momentum driven wind model (L 15N128csf,mdvwl .25) and according to t] ~ v„ 2 for the energy driven wind model 
(L15N128csf,edvwl.25). A fit to the observed GSMF llLi fe Whitdl200^ is shown for reference. Also shown are constraints on halo stellar 
fractions obtained using the abundance matching technique llGuo et alj feoiO). These results illustrate that the overproduction of low 
mass galaxies is reduced in models in which the mass loading increases for low escape velocities. 




Figure 3. Dependence of the low-mass end of the GSMF (left panel) in energy driven wind models (n ~ v^ 2 ) on the assumed normaliza- 
tion of wind velocity and mass loading. The effects on the stellar mass fractions of halos are illustrated in the figure's right panel. Shown 
are results for z = and models in which the wind velocity is given by v w = 1.25 X v esc (L15N128csf,edvwl.25) and v w = 0.6 X v esc 
(L15N128csf,edvw0.6). In the latter models, winds have an accordingly larger mass loading as n ~ v^, 2 . Also illustrated are the effects of 
increasing the available supernova feedback energy by assuming a Chabrier rather than a Salpeter IMF and of assuming more efficient 
photoheating (soli d green line, simulation L15N128csf,edvw0.6,uv,chab, see main text for more details). A fit to the observed GSMF 
dLi fc Whitell2009h is shown fo r reference. Also shown are constraints on halo stellar fractions obtained using the abundance match- 
ing technique l|Guo et al.ll2010h . These results illustrate that winds with lower velocity but larger mass loading are more successful in 
suppressing the overproduction of low mass galaxies. As expected, boosting supernova feedback and photoheating also reduces their 
number. 
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at z = 0. As a next step, we will explore whether this result 
is sensitive to changes in the simulation resolution. Figure [4] 
shows the GSMFs and stellar mass fractions for runs with- 
out winds, with constant winds, as well as for the aforemen- 
tioned energy-driven variable wind model. Results are given 
for our standard resolution simulations and high-resolution 
runs with eight times better mass resolution and a two times 
smaller gravitational softening. For the energy-driven wind 
model, we also performed a very high-resolution run with 
64 times better mass resolution and a four times smaller 
gravitational softening length. 

Without winds, the results are well converged with 
respect to resolution above a limiting mass of ~ 5 x 
10 9 h~ 2 M Q . Below this threshold, the GSMF in the high- 
resolution run keeps rising for decreasing stellar mass, while 
for the standard resolution it reaches a maximum and then 
falls off. This is a typical behaviour found in resolution stud- 
ies of cosmic structure formation simulations. The runs with 
constant winds exhibit a very similar behaviour, albeit at 
lower stellar mass values. 

Interestingly, the resolution limit in terms of halo mass 
is very similar in runs without and with constant winds, 
roughly lying at ~ 7 x 1O 1O M0, while it is, due to the much 
smaller stellar fractions, shifted towards lower stellar mass 
values in the latter runs. This suggests that accurately re- 
solving the host halo and its potential well is more important 
than the number of star particles which represent a galaxy, 
at least for the quantities analysed here. 

For the variable wind simulations, the situation is some- 
what different. The slope of the low-mass end of the observed 
GSMF is surprisingly well reproduced at all three resolution 
levels even at very small stellar mass. The normalisation, 
however, increases somewhat between the standard and high 
resolution runs. Compared to the latter, the normalisation 
is only slightly larger in the very high resolution simula- 
tion, most likely indicating that we are approaching conver- 
gence. On second thought, it is not entirely unexpected that 
the convergence behaviour is different in the variable wind 
model. In low-mass galaxies it implies much smaller wind 
velocities than the constant wind scenario. Consequently, 
the maximum radius that is reached by an outflow is much 
more sensitive to the detailed structure of the galaxy, which 
in turn depends on the gravitational softening and the mass 
resolution in the simulation. 

Overall, the galaxy populations in the simulations with- 
out and with constant winds seem to be well converged 
above the quoted resolution limits. They however fail to 
match the observed low-mass end of the GSMF. In the 
energy-driven variable wind runs, the low-mass end slope 
of the GSMF is reproduced well, while the normalisation is 
close to converged in the two higher resolution runs in the 
stellar mass range shown in Fig. [4] 



3.2 The high-mass end of the galaxy stellar mass 
function 

The left panel of Fig.f3]suggests that supernova-driven winds 
do not strongly affect the abundance of large galaxies with 
stellar masses above ~ 10 11 /i _1 tf0. In the constant wind 
model the wind velocity becomes insufficient to drive gas 
to large radii or to expel it from these objects, while in 



the variable wind model the wind's mass loading becomes 
increasingly small for very massive galaxies. 

Feedback from AGN, instead, is widely considered to 
be a promising candidate for suppressing star formation 
in massive galaxies. In particular, it appears to be a nec- 
essary i ngredient for explaining the high-m ass end of the 
GSMF (|Croton et al.ll2006l ; [Bower et al.l2006f ). Another, po- 
tentially indispensable ingredient that has been ignored in 
many previous numerical studies is to properly account for 
diffuse stellar light. A stellar component with very low sur- 
face brightness might be missed in observations or not be 
considered to be part of a galaxy. However, when GSMFs 
are measured from simulations in a straightforward way the 
mass of these stars is normally added to the stellar mass of 
the galaxy that resides at the centre of the corresponding 
halo. Simply because of a different accounting of the diffuse 
light, this can induce a discrepancy between observed and 
simulated GSMFs, particularly at the high-mass end, where 
the fraction of diffuse stars can be quite large. 

Figure [5] illustrates how accounting for AGN feedback 
and diffuse stellar light affects the high-mass end of the 
GSMF. AGN feedback is included as described in Sec- 
tion 12.31 A correction for diffuse stellar light is included 
based on a surface brightness threshold, i.e. stars that are 
seen as diffuse light dimmer than 26.5 mag arcsec -2 in the 
V-band or 20.7 mag arcsec -2 in the K-band are not assigned 
to any galaxy and ignored when computing the GSMF. Fig- 
ure [S] clearly shows that both AGN feedback and the cor- 
rection for a diffuse stellar component reduce the overabun- 
dance of massive galaxies substantially, thereby helping to 
reproduce the high-mass end of the observed GSMF. 

The AGN model employed here efficiently suppresses 
late-time star form ation in massive galaxies (see also 
iPuchwein et al.l |2010T ). However, its exact impact on the 
GSMF depends somewhat on the adopted model parame- 
ters. For the choices outlined in Section T2.3I we get a quite 
good fit to the observed GSMF. We also did some tests with 
other parameter values and found that accurately reproduc- 
ing the knee of the GSMF is sensitive to the value of the 
seed BH mass. More specifically, choosing a larger value for 
the seed black hole mass while keeping the threshold mass of 
the FoF groups into which they are implanted unchanged re- 
sults in AGN that become too efficient in lower-mass galax- 
ies. The abundance of galaxies around the observed knee 
of the GSMF is then underpredicted. Conversely, lowering 
the seed BH mass is accordingly expected to cause AGN to 
become important only for more massive galaxies, thereby 
shifting the knee of the simulated GSMF to larger stellar 
mass values. Note that the the position of the knee also de- 
pends on the exact choice of the AGN feedback parameters 
in semi-analytical galaxy formation models. There, typically 
a threshold for the ratio of cooling to free-fall time is used 
to control essentially the halo mass above which AGN feed- 
back is effective, t hereby shifting the knee of the GSMF (e.g. 
iBower et af]|2012l ). 

Correcting for a dim diffuse stellar component lowers 
the predicted abundance of galaxies more and more strongly 
for increasing stellar mass. This trend can be understood 
as a result of the fact that central galaxies in groups and 
clusters have particularly large diffuse light envelopes. Also 
note that we do not exclude stars from poorly resolved low- 
mass galaxies that are within one softening length from the 
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Figure 4. Dependence of the low-mass end of the GSMF (left panel) and the stellar mass fractions of halos (right panel) at z = 
on the simulation resolution. The effects on the stellar mass fractions of halos are illustrated in the figure's right panel. Shown 
are results obtained without kinetic wind feedback (L15N128csf, L15N256csf), with winds with constant velocity and mass loading 
(L15N128csf,cw, L15N256csf,cw), and with energy driven variable winds (L15N128csf,edvw0.6,uv,chab, L15N256csf,edvw0.6,uv,chab, 
L15N512csf,edvw0.6,uv,chab). All simulations shown here adopt a box size of 15 h~ 1 Mpc. Results are given for simulations with 2 X 128 3 , 
2 X 256 3 , and 2 X 512 3 particles (solid, dashed, and dotted lines, respectively), where the latter resolution is only available for the energy- 
driven variable wind model. 
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Figure 5. The impact of AGN feedback and diffuse light subtraction on the high-mass end of the GSMF (left panel) and the stellar 
mass fractions of halos (right panel) at z = 0. The results are based on runs with energy driven variable winds and were performed with 
(red curves, L60N512csfbh,edvw0.6,uv,chab) and without black holes (green curves, L60N512csf,edvw0.6,uv,chab). We either accounted 
for all stars in a galaxy's host halo (solid lines) or subtracted a diffuse stellar component defined by having a V-band surface brightness 
dimmer than 26.5 magarcsec -2 (dotted) or a K-band surface brightness dimmer than 20.7 magarcsec -2 (dashed). Both, accounting for 
AGN feedback and excluding diffuse light, significant ly reduces the overproduction of massive galaxies and improves the agreement with 
observations. The observed GSMF jLi k Whitdl2009l . squares with error bars) and a Schechter fit to it are shown for referenc e in the left 
pane l, while the right panel displays constraints on halo stellar fractions obtained using the abundance matching technique l lGuo et alj 
l20irl for comparison. 



galaxy's centre. As expected, the exact change in the high- 
mass end of the GSMF depends on the adopted surface 
brightness threshold value. For an accurate comparison of 
simulations and observations it is thus important that ob- 
servational studies report the surface brightness level down 
to which stars are accounted for in observational constraints 



on the GSMF. Simulators then need to adopt this threshold 
in their analysis. 

To summarize, the high-mass end of the GSMF can 
be reproduced quite well when both AGN feedback and a 
correction for diffuse light is included. In combination with 
our energy-driven variable wind model, the predicted GSMF 
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Figure 6. Star formation rate density as a function of redshift 
for different simulation models. Results are shown for simula- 
tions without winds (L15N256csf), as well as for runs with con- 
stant winds (L15N256csf,cw) and energy-driven variable winds 
(L15N256csf,edvw0.6,uv,chab). For the latter wind model also 
predictions based on a run with AGN feedback are shown 
(L60N512csfbh,edvw0.6,uv,chab). The dark grey lines indicate 
different fits to observational data (some of the data is also indi- 
cated by squares with error bars, see main text for more details). 

then matches the observations impressively well over the 
whole stellar mass range shown in Fig. [5] 



3.3 The cosmic star formation history and high 
redshift stellar mass functions 

In the previous sections, we have identified a feedback model 
that reproduces the GSMF at redshift z = well. We will 
now explore how the cosmic star formation rate history and 
high-redshift GSMFs are affected by winds and AGN, and 
whether the predictions of this favoured feedback model are 
compatible with observations. 

Figure [6] shows the star formation rate density as a 
function of redshift for runs without winds, as well as for 
simulations with constant and energy-driven variable winds. 
For the latter model, we also include results for a run that 
additionally incl udes AGN f eedba ck. F its to observational 
data obtained in ICole et all (|200lD and|Hopk ins fc Beacoml 
(2006) are shown for reference. In the latter case, they are 
displayed both fo r a modified Salpeter and a Baldry & 
Glazebrook IMF l|Baldrv fc Glazebrookl 120031 ) . For the re- 
sults based on the Salpeter IMF also the data including their 
errors are indicated in the figure. The large error bars and 
significant scatter in the data also illustrate that the star 
formation rate density, especially at high redshift, is notori- 
ously difficult to measure, in particular due to uncertainties 
in the IMF and the dust correction. 

Energy-driven variable winds can efficiently suppress 
star formation at high redshift where it happens mostly in 
small objects, while runs without winds or with constant 
winds overpredict the cosmic star formation rate density for 
z > 1 and z > 4, respectively. In addition, the peak of the 



star formation rate density shifts to lower redshift (z ~ 2.5) 
in the variable wind model, in good agreement with obser- 
vations. In the other models, instead, the star formation his- 
tory peaks at higher redshift, which seems disfavoured by the 
data. As a side note, the small "wiggles" in the curves of the 
variable wind models correspond to the times at which the 
on-the-ffy FoF group finder is run and at which the host halo 
masses of all gas particles in the simulation, which are used 
to compute wind velocities and mass loadings, are updated. 

In most of the redshift range shown in Fig. [6] the 
star formation rate density in our energy-driven variable 
wind simulations is slightly lower than suggested by obser- 
vations. At the same time, the low-redshift GSMF seems to 
be slightly overpredicted. This suggests tensions between the 
observational constraints on the redshift evolution of the star 
formation rate density and the ob served low-redshift G SMF 
as previously pointed out by, e.g.. iRudnick et all (|2003l ) and 
iHopkins fc Beacoml (|2006l ). 

At low redshift the observed star formation rate density 
drops almost an order of magnitude below its peak value. A 
similarly strong decrease can be seen in runs without winds. 
There it is, however, a consequence of the almost complete 
consumption of gas in galaxies and is closely connected to ex- 
cessively large stellar mass fractions. In constant wind sim- 
ulations the star formation rate density falls off too early, 
and is smaller than observed for z < 3. The decline is how- 
ever less steep, so that at z = it becomes comparable to 
the observed values again. The star formation rate density 
in the variable wind model also decreases at low redshift. 
While it remains roughly comparable to the observed values 
at all times, the decline seems to be somewhat less steep 
than observed. 

We also investigated how AGN feedback affects the cos- 
mic star formation history and found that it does not have 
a huge impact at high redshift but results in a steep drop-off 
below z — 1. While a strong decline at low redshift is, in- 
deed, suggested by observations, this decline seems to start 
too early in our run with AGN. 

However, further test simulations showed that the onset 
of this drop-off is sensitive to the exact choice made for the 
AGN model parameters, similar to what we found for the 
position of the knee of the GSMF. In particular, choosing 
a larger seed BH mass shifts the decline to earlier times. 
Thus, we expect that it should also be possible to shift it to 
the observed later time by adopting a correspondingly tuned 
smaller value. 

Figure [7] displays the GSMFs at redshifts z = 1 and 2 
for runs without winds, with constant winds and with en- 
ergy driven variable winds. Fits to the observed GSMF from 
ISantini et all l|2012t ) are shown for comparison. These con- 
straints are slightly offset in redshift as they were not avail- 
able for the exact times of our simulation outputs. Clearly, 
the simulation without winds overpredicts the GSMF at 
z = 1 and 2. The constant wind model is in reasonably 
good agreement with the z — 2 GSMF over the limited stel- 
lar mass range in which it is constrained. At z = 1, however, 
the slope of the GSMF is much steeper than observed. On 
the other hand, the energy-driven variable wind model is in 
very good agreement with the data both at redshifts z = 1 
and 2. 

Overall, the variable wind model, whose parameter 
choice was guided by data at z ~ 0, also matches the 
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Figure 7. The impact of galactic winds on the z = 1 (left panel) and z = 2 (right panel) GSMF. Results are shown for runs without 
winds (L60N512csf), with c onstant winds (L60N 512csf,cw), and with energy driven variable winds (L60N512csf,edvw0.6,uv,chab). Fits 
to the observed GSMF from lSantini et al.l j2012D are shown for comparison. The black solid lines indicate the stellar mass range in which 
the fits are constrained by the data, while the black dotted lines are extrapolations. Also at z = 1 and 2, simulations with energy-driven 
variable winds are in very good agreement with the data. 



observed cosmic star formation history and high-redshift 
GSMFs very well. In contrast, the constant and no wind 
models clearly fail in doing so. 



Specific star formation rates and gas fractions 
of galaxies 



3.4 



Our analysis in Sections 13. II and !3. 21 was based on the stellar 
masses of galaxies at z = 0, and thus essentially probed the 
integral of the star formation rate form early times to the 
present day. In the following we will, instead, examine the 
instantaneous star formation rates of galaxies at z — and 
their dependence on stellar mass. We will also study the gas 
fractions of galaxies which are of course closely related to 
their star formation rates. 

Figure [5] shows the specific star formation rate, i.e. the 
star formation rate per unit stellar mass, as a function of 
the stellar mass of the galaxy at z — 0. Fits to observational 
data obtained bv lSalim et al.l (|2007h are shown for reference. 
Clearly, the constant and no wind models fail to account 
for the observed specific star formation rates of galaxies, 
while the energy-driven variable wind model is in impres- 
sively good agreement with the data. 

This can be understood by looking at the amount of 
cold gas in galaxies, or in other words the available fuel for 
star formation. Figure[JJ]displays the gas mass to stellar mass 
ratio of galaxies as a function of their stellar mass. As we 
are interested in cold gas rather than in hot halo gas, and in 
order to make our results directly comparable to Hi obser- 
vations, we included only neutral hydrogen when computing 
the gas masses of simulated galaxies. To this end we assumed 
that the star-forming multiphase interstellar medium with 
hydrogen densities nn > 0.23 cm -3 is neutral. In the em- 
ployed star formation model, this is a reasonable assumption 
as, somewhat depending on the exact density, at least 84% 
of the ISM are supposed to be in a cold phase with tempera- 
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Figure 8. Specific star formation rates (SSFRs), i.e. star forma- 
tion rates per unit stellar mass, for galaxies simulated without 
winds (L15N256csf), with constant winds (L15N256csf,cw), and 
with energy driven variable winds (L15N 256csf,edvw0 . 6,uv ,chab) 
at z = 0. Fits to observed SSFRs from I Salim et al] 1 120071 ) are 
displayed for reference. The solid and dashed grey lines indi- 
cate measurements for "pure" star-forming galaxies and star- 
forming galaxy/AGN composites, respectively. The energy-driven 
wind model produces galaxies with much larger SSFRs, thereby 
improving the agreement with observational constraints signifi- 
cantly. 



tures below T = 10 3 K (see lSpringel fc Hernau"isl2003r ). For 
the remaining cold gas that falls below this density threshold 
and that is part of a galaxy's Subfind (sub-)halo, the neu- 
tral hydrogen fraction is computed explicitly as a function 
of gas temperature and density under the assumption of col- 
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Figure 9. Mean gas mass (neutral hydrogen) to stellar mass 
ratios of galaxies simulated without winds (L15N256csf ), with 
constant winds (L15N256csf,cw), and with energy driven vari- 
able winds (L15N256csf,edvw0.6,uv,chab), as a function of stel- 
lar mass at z = 0. Bi nned observational data compiled by 
iPeeples fe Shankarl ll201ll) are shown for comparison. The error 
bars indicate the standard deviations of the bins' mean values. 
The energy-driven variable wind model produces galaxies with 
much larger gas mass to stellar mass ratios, bringing them into 
excellent agreement with observational constraints. 

lisional ionisation equi librium in the pre sence of an external 
UV background, as in lKatz et al.l (| 19961 ). 

The gas mass to stellar mass ratios exhibit a similar 
behaviour as the specific star formation rates. Runs with- 
out and with constant winds underpredict the observed gas 
mass to stellar mass ratios, while the predictions based on 
the energy-driven variable wind model are in excellent agree- 
ment with the data. This is independent evidence that the 
latter model ejects just the right amount of gas from galax- 
ies, enabling them to simultaneously reproduce both the 
observed stellar and gas masses. The gas fractions of the 
constant and no wind models look surprisingly similar. The 
cause for the low gas content is, however, different in these 
models. In the constant wind run, outflows expel a large 
fraction of the gas from galaxies, except in very massive ob- 
jects with stellar masses above ~ 10 11 Mq, in which the fixed 
wind velocity becomes insufficient to unbind gas. Without 
winds, the gas fractions are instead reduced by a highly ef- 
ficient conversion of gas into stars. It thus appears to be a 
coincidence that the curves are so similar below ~ 1O 11 M0 
at z = 0. This is also suggested by Fig. [6] which shows that 
the star formation rate history is quite different in these 
runs, only at 2 = the star formation rate density happens 
to be remarkably similar. 

These results demonstrate that energy-driven variable 
winds bring not only the stellar masses and star formation 
histories of galaxies into much better agreement with obser- 
vations but also yield gas masses and current specific star 
formation rates in good agreement with observational con- 
straints. We note that this impressively good agreement with 
the data is found over the whole stellar mass range consid- 
ered here. 



4 SUMMARY AND CONCLUSIONS 

In this work, we investigated feedback processes relevant for 
shaping the GSMF, and how they can be accounted for in 
cosmological hydrodynamical simulations of galaxy forma- 
tion as a subresolution model. For this purpose, we imple- 
mented a new model for galactic winds into the Gadget- 
3 simulation code and used it to perform a set of high- 
resolution simulations. This allowed us to study how wind 
velocity and wind mass flux need to scale with the proper- 
ties of the galaxy from which the winds are launched such 
that the resulting simulated galaxy populations match ob- 
servations. Even though identifying a favoured model in this 
fashion is phenomenological in nature, we emphasize that a 
successful wind model is a great asset for future galaxy sim- 
ulations and attempts to understand the detailed physical 
processes that launch winds with the required properties in 
the first place. 

We also studied the impact of feedback from AGN on 
the abundance of massive galaxies and assessed the impor- 
tance of correctly accounting for diffuse light when compar- 
ing this quantity between simulations and observations. Our 
main findings are: 

• Galactic winds with constant mass-loading and velocity, 
which have been widely used in the literature, are unable to 
reproduce the observed shape and shallow faint-end slope of 
the GSMF. 

• Adopting a variable wind velocity equal to a specific 
fraction of the escape velocity of the galaxy, as suggested 
by observations, and allowing either for a certain amount of 
momentum or energy to drive the wind, much larger wind 
mass-loadings in low-mass galaxies are possible. This results 
in a shallower low-mass end of the GSMF. 

• A stronger scaling of mass-loading with wind velocity, 
as implied by an energy- compared to a momentum-driven 
wind scenario, reduces the slope of the low-mass end of the 
GSMF more strongly. We find that an energy-driven variable 
velocity wind model reproduces the observed low-mass slope 
well. 

• In such an energy-driven model the normalisation of the 
GSMF depends on the available energy for driving the wind, 
as well as on the adopted ratio of wind velocity to escape 
velocity. 

• The observed normalisation of the low-mass-end of 
the GSMF can be matched well when wind velocities are 
adopted that are somewhat smaller than the nominal escape 
velocity from the centre of the halo (in agreement with ob- 
servational suggestions) , and when calculating the available 
energy for driving the wind based on the expected supernova 
rate. 

• Including AGN feedback and correcting for diffuse light 
when computing galaxy masses reduces the abundance of 
massive galaxies, thereby bringing the high-mass end of the 
GSMF into much better agreement with observations. 

Overall, our simulations with energy-driven variable 
winds and AGN feedback reproduce the observed GSMF 
at redshift z = quite well. In order to gain further insights 
in how viable such scenarios are in light of other data, we 
investigated how simulations with these models compare to 
observational constraints on the cosmic star formation his- 
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tory, as well as with observations of the GSMFs at redshifts 
z — 1 and 2. We found that: 

• Energy-driven variable winds suppress star formation 
at high redshift, reducing the cosmic star formation rate 
density to observed levels and shifting its peak to z ~ 2.5, 
which is in good agreement with the data. 

• The cosmic star formation rate density declines at low 
redshift in the energy-driven wind model. The drop-off, how- 
ever, seems to be somewhat less steep than observed, but 
steepens when AGN feedback is also included. 

• When the whole redshift range of < z < 5 is consid- 
ered, the energy-driven variable wind model reproduces the 
overall shape and normalisation of the observed cosmic star 
formation rate history much better than any of the other 
models. 

• The observed GSMFs at redshifts z = 1 and 2 are repro- 
duced very well by the energy-driven variable wind model, 
whereas simulations without winds fail to do so. Especially 
at z = 1, runs with constant mass- loading and wind velocity 
also fail. 

• Over the whole stellar mass range considered, we find 
that the cold gas masses and specific star formation rates at 
z — are in excellent agreement with observations in runs 
with energy-driven variable winds. 

• In contrast, in simulations without winds or with fixed 
wind velocity and mass loading, the gas mass to stellar mass 
ratios and specific star formation rates are underpredicted 
by almost one order of magnitude. 

To summarise, our new energy-driven variable velocity 
wind model faithfully reproduces the observed low-mass end 
of the galaxy stellar mass function. For the first time, wc 
combined such variable winds with a model for AGN feed- 
back and with a proper treatment of diffuse light. Account- 
ing for these effects allows us to simultaneously reproduce 
the low-mass end, the knee region, and the high-mass end of 
the observed GSMF. We also obtained stellar mass fractions 
of halos that compare very well with results from abundance 
matching. 

Furthermore, this model does not only make the stellar 
masses of galaxies compatible with observations but simul- 
taneously almost perfectly matches their gas mass to stellar 
mass ratios as well as their specific star formation rates. In 
addition, the model drastically improves the agreement with 
observed GSMFs at high redshift, and with constraints on 
the cosmic star formation history. 

This broad agreement with a wide range of data is a sig- 
nificant improvement compared to earlier attempts to model 
the buildup of the whole galaxy population in our Universe 
through cosmological hydrodynamical simulations. For the 
first time, such hydrodynamic simulations achieve a com- 
parable success as semi-analytic models. Simulation models 
that reproduce the main properties of the observed galaxy 
population in this fashion will be extremely helpful for mak- 
ing more realistic simulation predictions in the future, for 
example for the clustering of galaxies, for the galaxy bias, 
for chemical enrichment, for galaxy merger rates as a func- 
tion of stellar mass, for the formation of intracluster light, 
etc., to just name a few. 

Also, while the present work provides no insights into 
the detailed physics of the creation of winds in galaxies, our 
results strongly suggest that the mass loading in real galactic 



winds needs to strongly increase in low-mass galaxies, as 
otherwise an overproduction of the faint-end of the galaxy 
luminosity function seems inevitable. At the same time, our 
work demonstrates that the energy available from feedback 
processes associated with star formation is sufficient to drive 
such winds if losses by radiative cooling are not dominant. 
These findings can provide some guidance in the ongoing 
quest to fully understand wind formation in galaxies. 
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